Comparative Transcriptome Analysis of Babesia bigemina Attenuated Vaccine and Virulent Strains of Mexican Origin.

Bovine babesiosis, caused by the protozoan Babesia bigemina, is one of the most important hemoparasite diseases of cattle in Mexico and the world. An attenuated B. bigemina strain maintained under in vitro culture conditions has been used as a live attenuated vaccine; however, the biological mechanisms involved in attenuation are unknown. The objective of this study was to identify, through a comparative transcriptomics approach, the components of the B. bigemina virulent parasites that are differentially expressed in vivo, as opposed to those expressed by B. bigemina attenuated vaccine parasites when inoculated into naïve cattle. The biological material under study was obtained by inoculating spleen-intact cattle with infected erythrocytes containing either the attenuated strain or a virulent field strain. After RNA extraction, transcriptomic analysis (RNA-seq) was performed, followed by bioinformatic Differential Expression (DE) analysis and Gene Ontology (GO) term enrichment. The high-throughput sequencing results obtained by analyzing three biological replicates for each parasite strain ranged from 9,504,000 to 9,656,000, and 13,400,000 to 15,750,000 reads for the B. bigemina attenuated and virulent strains, respectively. At least 519 differentially expressed genes were identified in the analyzed strains. In addition, GO analysis revealed both similarities and differences across the three categories: cellular components, biological processes, and molecular functions. The attenuated strain of B. bigemina derived from in vitro culture presents global transcriptomic changes when compared to the virulent strain. Moreover, the obtained data provide insights into the potential molecular mechanisms associated with the attenuation or pathogenicity of each analyzed strain, offering molecular markers that might be associated with virulence or potential vaccine candidates.

Exploration of low-frequency allelic variants of SARS-CoV-2 genomes reveals coinfections in Mexico occurred during periods of VOCs turnover.

A total of 14 973 alleles in 29 661 sequenced samples collected between March 2021 and January 2023 by the Mexican Consortium for Genomic Surveillance (CoViGen-Mex) and collaborators were used to construct a thorough map of mutations of the Mexican SARS-CoV-2 genomic landscape containing Intra-Patient Minor Allelic Variants (IPMAVs), which are low-frequency alleles not ordinarily present in a genomic consensus sequence. This additional information proved critical in identifying putative coinfecting variants included alongside the most common variants, B.1.1.222, B.1.1.519, and variants of concern (VOCs) Alpha, Gamma, Delta, and Omicron. A total of 379 coinfection events were recorded in the dataset (a rate of 1.28 %), resulting in the first such catalogue in Mexico. The most common putative coinfections occurred during the spread of Delta or after the introduction of Omicron BA.2 and its descendants. Coinfections occurred constantly during periods of variant turnover when more than one variant shared the same niche and high infection rate was observed, which was dependent on the local variants and time. Coinfections might occur at a higher frequency than customarily reported, but they are often ignored as only the consensus sequence is reported for lineage identification.

Multiomics analysis reveals B. MO1 as a distinct Babesia species and provides insights into its evolution and virulence.

Babesiosis, caused by protozoan parasites of the genus Babesia , is an emerging tick-borne disease of significance for both human and animal health. Babesia parasites infect erythrocytes of vertebrate hosts where they develop and multiply rapidly to cause the pathological symptoms associated with the disease. The identification of various Babesia species underscores the ongoing risk of new zoonotic pathogens capable of infecting humans, a concern amplified by anthropogenic activities and environmental shifts impacting the distribution and transmission dynamics of parasites, their vectors, and reservoir hosts. One such species, Babesia MO1, previously implicated in severe cases of human babesiosis in the midwestern United States, was initially considered closely related to B. divergens , the predominant agent of human babesiosis in Europe. Yet, uncertainties persist regarding whether these pathogens represent distinct variants of the same species or are entirely separate species. We show that although both B. MO1 and B. divergens share similar genome sizes, comprising three nuclear chromosomes, one linear mitochondrial chromosome, and one circular apicoplast chromosome, major differences exist in terms of genomic sequence divergence, gene functions, transcription profiles, replication rates and susceptibility to antiparasitic drugs. Furthermore, both pathogens have evolved distinct classes of multigene families, crucial for their pathogenicity and adaptation to specific mammalian hosts. Leveraging genomic information for B. MO1, B. divergens , and other members of the Babesiidae family within Apicomplexa provides valuable insights into the evolution, diversity, and virulence of these parasites. This knowledge serves as a critical tool in preemptively addressing the emergence and rapid transmission of more virulent strains.

SARS-CoV-2 Omicron variants BA.4 and BA.5 dominated the fifth COVID-19 epidemiological wave in Mexico.

In Mexico, the BA.4 and BA.5 Omicron variants dominated the fifth epidemic wave (summer 2022), superseding BA.2, which had circulated during the inter-wave period. The present study uses genome sequencing and statistical and phylogenetic analyses to examine these variants' abundance, distribution, and genetic diversity in Mexico from April to August 2022. Over 35 % of the sequenced genomes in this period corresponded to the BA.2 variant, 8 % to the BA.4 and 56 % to the BA.5 variant. Multiple subvariants were identified, but the most abundant, BA.2.9, BA.2.12.1, BA.5.1, BA.5.2, BA.5.2.1 and BA.4.1, circulated across the entire country, not forming geographical clusters. Contrastingly, other subvariants exhibited a geographically restricted distribution, most notably in the Southeast region, which showed a distinct subvariant dynamic. This study supports previous results showing that this region may be a significant entry point and contributed to introducing and evolving novel variants in Mexico. Furthermore, a differential distribution was observed for certain subvariants among specific States through time, which may have contributed to the overall increased diversity observed during this wave compared to the previous ones. This study highlights the importance of sustaining genomic surveillance to identify novel variants that may impact public health.

Draft genome sequences of six bacteria isolated from polylactic acid samples from the Mediterranean Sea.

Here, we report the draft genome sequences of six marine strains isolated from plastic samples incubated in the Mediterranean Sea. Genomic analyses place these strains within the Alkalihalobacillus, Bacillus, Halomonas, and Marinobacter genera. Examining the genomes of these non-typical environmental bacteria increases our comprehension of microorganism biology and their potential uses.

Identification of Cry toxin receptor genes homologs in a de novo transcriptome of Premnotrypes vorax (Coleoptera: Curculionidae).

The white potato worm Premnotrypes vorax (Hustache) (Coleoptera: Curculionidae) is one of the most destructive insect pests of potato crops in South America. Like many coleopteran insects, P. vorax shows low susceptibility to Cry insecticidal proteins produced by the bacterium Bacillus thuringiensis (Bt). However, the presence of Cry toxin receptors in the midgut of this this insect has never been studied. The main Cry-binding proteins described in other insect species are cadherin (CAD), aminopeptidase N (APN), alkaline phosphatase (ALP) and ATP-binding cassette (ABC) transporters. In this study, we analyzed and validated a de novo assembled transcriptome of Illumina sequencing data to identify and to characterize homologs of Cry toxin receptors. We identified the protein sequences in P. vorax that show high identity with their orthologous sequences of the Cry toxin binding proteins in other coleopteran larvae such as APN, ALP, CAD and ABC transporter. This study provides preliminary identification of putative receptor genes of Cry proteins that would be useful for future studies involving biocontrol of this important potato crop pest.

Multiple approaches of cellular metabolism define the bacterial ancestry of mitochondria.

We breathe at the molecular level when mitochondria in our cells consume oxygen to extract energy from nutrients. Mitochondria are characteristic cellular organelles that derive from aerobic bacteria and carry out oxidative phosphorylation and other key metabolic pathways in eukaryotic cells. The precise bacterial origin of mitochondria and, consequently, the ancestry of the aerobic metabolism of our cells remain controversial despite the vast genomic information that is now available. Here, we use multiple approaches to define the most likely living relatives of the ancestral bacteria from which mitochondria originated. These bacteria live in marine environments and exhibit the highest frequency of aerobic traits and genes for the metabolism of fundamental lipids that are present in the membranes of eukaryotes, sphingolipids, and cardiolipin.

Comparing the evolutionary dynamics of predominant SARS-CoV-2 virus lineages co-circulating in Mexico.

Over 200 different SARS-CoV-2 lineages have been observed in Mexico by November 2021. To investigate lineage replacement dynamics, we applied a phylodynamic approach and explored the evolutionary trajectories of five dominant lineages that circulated during the first year of local transmission. For most lineages, peaks in sampling frequencies coincided with different epidemiological waves of infection in Mexico. Lineages B.1.1.222 and B.1.1.519 exhibited similar dynamics, constituting clades that likely originated in Mexico and persisted for >12 months. Lineages B.1.1.7, P.1 and B.1.617.2 also displayed similar dynamics, characterized by multiple introduction events leading to a few successful extended local transmission chains that persisted for several months. For the largest B.1.617.2 clades, we further explored viral lineage movements across Mexico. Many clades were located within the south region of the country, suggesting that this area played a key role in the spread of SARS-CoV-2 in Mexico.

Babesia duncani multi-omics identifies virulence factors and drug targets.

Babesiosis is a malaria-like disease in humans and animals that is caused by Babesia species, which are tick-transmitted apicomplexan pathogens. Babesia duncani causes severe to lethal infection in humans, but despite the risk that this parasite poses as an emerging pathogen, little is known about its biology, metabolic requirements or pathogenesis. Unlike other apicomplexan parasites that infect red blood cells, B. duncani can be continuously cultured in vitro in human erythrocytes and can infect mice resulting in fulminant babesiosis and death. We report comprehensive, detailed molecular, genomic, transcriptomic and epigenetic analyses to gain insights into the biology of B. duncani. We completed the assembly, 3D structure and annotation of its nuclear genome, and analysed its transcriptomic and epigenetics profiles during its asexual life cycle stages in human erythrocytes. We used RNA-seq data to produce an atlas of parasite metabolism during its intraerythrocytic life cycle. Characterization of the B. duncani genome, epigenome and transcriptome identified classes of candidate virulence factors, antigens for diagnosis of active infection and several attractive drug targets. Furthermore, metabolic reconstitutions from genome annotation and in vitro efficacy studies identified antifolates, pyrimethamine and WR-99210 as potent inhibitors of B. duncani to establish a pipeline of small molecules that could be developed as effective therapies for the treatment of human babesiosis.

SARS-CoV-2 BW lineage, a fast-growing Omicron variant from southeast Mexico bearing relevant escape mutations.

PURPOSE: The swift expansion of the BW.1 SARS-CoV-2 variant coincided with a rapid increase of COVID-19 cases occurring in Southeast Mexico in October, 2022, which marked the start of Mexico's sixth epidemiological wave. In Yucatan, up to 92% (58 of 73) of weekly sequenced genomes between epidemiological week 42 and 47 were identified as either BW.1 or its descendant, BW.1.1 in the region, during the last trimester of 2022. In the current study, a comprehensive genomic comparison was carried out to characterize the evolutionary history of the BW lineage, identifying its origins and its most important mutations. METHODS: An alignment of all the genomes of the BW lineage and its parental BA.5.6.2 variant was carried out to identify their mutations. A phylogenetic and ancestral sequence reconstruction analysis with geographical inference, as well as a longitudinal analysis of point mutations, were performed to trace back their origin and contrast them with key RBD mutations in variant BQ.1, one of the fastest-growing lineages to date. RESULTS: Our ancestral reconstruction analysis portrayed Mexico as the most probable origin of the BW.1 and BW.1.1 variants. Two synonymous substitutions, T7666C and C14599T, support their Mexican origin, whereas other two mutations are specific to BW.1: S:N460K and ORF1a:V627I. Two additional substitutions and a deletion are found in its descending subvariant, BW.1.1. Mutations found in the receptor binding domain, S:K444T, S:L452R, S:N460K, and S:F486V in BW.1 have been reported to be relevant for immune escape and are also key mutations in the BQ.1 lineage. CONCLUSIONS: BW.1 appears to have arisen in the Yucatan Peninsula in Southeast Mexico sometime around July 2022 during the fifth COVID-19 wave. Its rapid growth may be in part explained by the relevant escape mutations also found in BQ.1.

Omicron-BA.1 Dispersion Rates in Mexico Varied According to the Regional Epidemic Patterns and the Diversity of Local Delta Subvariants.

PURPOSE: The Omicron subvariant BA.1 of SARS-CoV-2 was first detected in November 2021 and quickly spread worldwide, displacing the Delta variant. In this work, a characterization of the spread of this variant in Mexico is presented. METHODS: The time to fixation of BA.1, the diversity of Delta sublineages, the population density, and the level of virus circulation during the inter-wave interval were determined to analyze differences in BA.1 spread. RESULTS: BA.1 began spreading during the first week of December 2021 and became dominant in the next three weeks, causing the fourth COVID-19 epidemiological surge in Mexico. Unlike previous variants, BA.1 did not exhibit a geographically distinct circulation pattern. However, a regional difference in the speed of the replacement of the Delta variant was observed. CONCLUSIONS: Viral diversity and the relative abundance of the virus in a particular area around the time of the introduction of a new lineage seem to have influenced the spread dynamics, in addition to population density. Nonetheless, if there is a significant difference in the fitness of the variants, or if the time allowed for the competition is sufficiently long, it seems the fitter virus will eventually become dominant, as observed in the eventual dominance of the BA.1.x variant in Mexico.

Argonaute5 and its associated small RNAs modulate the transcriptional response during the rhizobia-Phaseolus vulgaris symbiosis.

Both plant- and rhizobia-derived small RNAs play an essential role in regulating the root nodule symbiosis in legumes. Small RNAs, in association with Argonaute proteins, tune the expression of genes participating in nodule development and rhizobial infection. However, the role of Argonaute proteins in this symbiosis has been overlooked. In this study, we provide transcriptional evidence showing that Argonaute5 (AGO5) is a determinant genetic component in the root nodule symbiosis in Phaseolus vulgaris. A spatio-temporal transcriptional analysis revealed that the promoter of PvAGO5 is active in lateral root primordia, root hairs from rhizobia-inoculated roots, nodule primordia, and mature nodules. Transcriptional analysis by RNA sequencing revealed that gene silencing of PvAGO5 affected the expression of genes involved in the biosynthesis of the cell wall and phytohormones participating in the rhizobial infection process and nodule development. PvAGO5 immunoprecipitation coupled to small RNA sequencing revealed the small RNAs bound to PvAGO5 during the root nodule symbiosis. Identification of small RNAs associated to PvAGO5 revealed miRNAs previously known to participate in this symbiotic process, further supporting a role for AGO5 in this process. Overall, the data presented shed light on the roles that PvAGO5 plays during the root nodule symbiosis in P. vulgaris.

Metagenomic analysis reveals differences in the co-occurrence and abundance of viral species in SARS-CoV-2 patients with different severity of disease.

BACKGROUND: SARS-CoV-2 infections have a wide spectrum of clinical manifestations whose causes are not completely understood. Some human conditions predispose to severe outcome, like old age or the presence of comorbidities, but many other facets, including coinfections with other viruses, remain poorly characterized. METHODS: In this study, the eukaryotic fraction of the respiratory virome of 120 COVID-19 patients was characterized through whole metagenomic sequencing. RESULTS: Genetic material from respiratory viruses was detected in 25% of all samples, whereas human viruses other than SARS-CoV-2 were found in 80% of them. Samples from hospitalized and deceased patients presented a higher prevalence of different viruses when compared to ambulatory individuals. Small circular DNA viruses from the Anneloviridae (Torque teno midi virus 8, TTV-like mini virus 19 and 26) and Cycloviridae families (Human associated cyclovirus 10), Human betaherpesvirus 6, were found to be significantly more abundant in samples from deceased and hospitalized patients compared to samples from ambulatory individuals. Similarly, Rotavirus A, Measles morbillivirus and Alphapapilomavirus 10 were significantly more prevalent in deceased patients compared to hospitalized and ambulatory individuals. CONCLUSIONS: Results show the suitability of using metagenomics to characterize a broader peripheric virological landscape of the eukaryotic virome in SARS-CoV-2 infected patients with distinct disease outcomes. Identified prevalent viruses in hospitalized and deceased patients may prove important for the targeted exploration of coinfections that may impact prognosis.

Superspreading event of Covid-19 in adolescents: is there a difference between the vaccinated and the unvaccinated?

Not available.

The New SARS-CoV-2 Variants and Their Epidemiological Impact in Mexico.

The COVID-19 disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus started its deadly journey into a global pandemic in Wuhan, China, in December 2019, where it was first isolated. Subsequently, multiple variants of the virus have been identified worldwide. In this review, we discuss the overall landscape of the pandemic in Mexico, including its most prevalent variants, their surveillance at a genomic level, and how they impacted the epidemiology of the disease. We also evaluate the heterologous vaccination in Mexico and how it may have influenced group immunity and helped mitigate the pandemic. Finally, we present an integrated view that could help us to understand the pandemic and serve as an example of the situation in Latin America.

Dominance of Three Sublineages of the SARS-CoV-2 Delta Variant in Mexico.

In this study, we analyzed the sequences of SARS-CoV-2 isolates of the Delta variant in Mexico, which has completely replaced other previously circulating variants in the country due to its transmission advantage. Among all the Delta sublineages that were detected, 81.5 % were classified as AY.20, AY.26, and AY.100. According to publicly available data, these only reached a world prevalence of less than 1%, suggesting a possible Mexican origin. The signature mutations of these sublineages are described herein, and phylogenetic analyses and haplotype networks are used to track their spread across the country. Other frequently detected sublineages include AY.3, AY.62, AY.103, and AY.113. Over time, the main sublineages showed different geographical distributions, with AY.20 predominant in Central Mexico, AY.26 in the North, and AY.100 in the Northwest and South/Southeast. This work describes the circulation, from May to November 2021, of the primary sublineages of the Delta variant associated with the third wave of the COVID-19 pandemic in Mexico and highlights the importance of SARS-CoV-2 genomic surveillance for the timely identification of emerging variants that may impact public health.

The Genomes of Two Strains of Taenia crassiceps the Animal Model for the Study of Human Cysticercosis.

Human cysticercosis by Taenia solium is the major cause of neurological illness in countries of Africa, Southeast Asia, and the Americas. Publication of four cestode genomes (T. solium, Echinococcus multilocularis, E. granulosus and Hymenolepis microstoma) in the last decade, marked the advent of novel approaches on the study of the host-parasite molecular crosstalk for cestode parasites of importance for human and animal health. Taenia crassiceps is another cestode parasite, closely related to T. solium, which has been used in numerous studies as an animal model for human cysticercosis. Therefore, characterization of the T. crassiceps genome will also contribute to the understanding of the human infection. Here, we report the genome of T. crassiceps WFU strain, reconstructed to a noncontiguous finished resolution and performed a genomic and differential expression comparison analysis against ORF strain. Both strain genomes were sequenced using Oxford Nanopore (MinION) and Illumina technologies, achieving high quality assemblies of about 107 Mb for both strains. Dotplot comparison between WFU and ORF demonstrated that both genomes were extremely similar. Additionally, karyotyping results for both strains failed to demonstrate a difference in chromosome composition. Therefore, our results strongly support the concept that the absence of scolex in the ORF strain of T. crassiceps was not the result of a chromosomal loss as proposed elsewhere. Instead, it appears to be the result of subtle and extensive differences in the regulation of gene expression. Analysis of variants between the two strains identified 2,487 sites with changes distributed in 31 of 65 scaffolds. The differential expression analysis revealed that genes related to development and morphogenesis in the ORF strain might be involved in the lack of scolex formation.

Association of Gut Microbiota with Dietary-dependent Childhood Obesity.

AIM: To evaluate the taxonomic profile of the gut microbiota using metagenomics and the association with diet-dependent childhood obesity. METHODS: A cross-sectional study of a subsample of 46 children was conducted. The children were classified as normal-weight, overweight, and obese according to their age and sex and the World Health Organization (WHO) guidelines. Dietary patterns were determined through principal component analysis. The profile of the human gut microbiota was determined by bioinformatic analysis using whole metagenome shotgun sequencing. The association of gut microbiota and z-BMI, waist circumference and hip circumference, and the possible modifying effect of diet were analyzed using multiple regression models. RESULTS: Children with an abundance of Holdemania spp. and high protein and complex carbohydrate consumption had a lower z-BMI (ß -19.06, p = 0.011), waist circumference (ß -171.92, p = 0.003), and hip circumference (ß -157.57, p = 0.004). In contrast, observed a positive association between Coprococcus catus and the low intake of this dietary pattern with hip circumference (ß 147.87, p = 0.025). Furthermore, the presence of Bilophila spp. and Paraprevotella xylaniphila with high saturated fat and simple carbohydrate consumption we observed a positive association between z-BMI (ß 47.5, p = 0.002), hip circumference (ß 44.54, p = 0.025), and waist circumference (ß 44.34, p = 0.004). CONCLUSION: We suggest that the synergism between diet and the profile of children's gut microbiota can be a factor that could be associated with the development of obesity and its complications in childhood.

The Alpha Variant (B.1.1.7) of SARS-CoV-2 Failed to Become Dominant in Mexico.

During the coronavirus disease 2019 (COVID-19) pandemic, the emergence and rapid increase of the B.1.1.7 (Alpha) lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), first identified in the United Kingdom in September 2020, was well documented in different areas of the world and became a global public health concern because of its increased transmissibility. The B.1.1.7 lineage was first detected in Mexico during December 2020, showing a slow progressive increase in its circulation frequency, which reached its maximum in May 2021 but never became predominant. In this work, we analyzed the patterns of diversity and distribution of this lineage in Mexico using phylogenetic and haplotype network analyses. Despite the reported increase in transmissibility of the B.1.1.7 lineage, in most Mexican states, it did not displace cocirculating lineages, such as B.1.1.519, which dominated the country from February to May 2021. Our results show that the states with the highest prevalence of B.1.1.7 were those at the Mexico-U.S. border. An apparent pattern of dispersion of this lineage from the northern states of Mexico toward the center or the southeast was observed in the largest transmission chains, indicating possible independent introduction events from the United States. However, other entry points cannot be excluded, as shown by multiple introduction events. Local transmission led to a few successful haplotypes with a localized distribution and specific mutations indicating sustained community transmission. IMPORTANCE The emergence and rapid increase of the B.1.1.7 (Alpha) lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) throughout the world were due to its increased transmissibility. However, it did not displace cocirculating lineages in most of Mexico, particularly B.1.1.519, which dominated the country from February to May 2021. In this work, we analyzed the distribution of B.1.1.7 in Mexico using phylogenetic and haplotype network analyses. Our results show that the states with the highest prevalence of B.1.1.7 (around 30%) were those at the Mexico-U.S. border, which also exhibited the highest lineage diversity, indicating possible introduction events from the United States. Also, several haplotypes were identified with a localized distribution and specific mutations, indicating that sustained community transmission occurred in the country.

Intra-host genetic population diversity: Role in emergence and persistence of drug resistance among Mycobacterium tuberculosis complex minor variants.

Drug resistant tuberculosis (DR-TB) is an important public health issue in different parts of the world. Mycobacterium tuberculosis complex variants (MTBC vars) preferentially infect certain hosts, limiting their distribution to different ecosystems. However, MTBC vars can infect other hosts beyond their preferred target potentially contributing to persistence of drug resistance (DR) in other niches. Here, we performed a comprehensive intra-host genetic analysis for the identification of DR-related mutations among all MTBC minor vars whole genome sequences (8,095 strains) publicly available worldwide. High confidence drug-resistance mutations in katG (isoniazid), rpsL (streptomycin), pncA (pyrazinamide), rpoB (rifampicin) and gyrA (fluoroquinolones) genes were identified among intrahost minor sub-populations in 197 different strains (2.43%) belonging to vars africanum, bovis, caprae, microti, orygis and pinnipedii. In addition, a three-dimensional structure modeling analysis to assess the role of novel mutations was also performed. Our findings highlight the importance of detecting discrete intra-host populations carrying DR mutations.